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ABSTRACT 

We explore the properties of dust and associated molecular gas in 352 nearby 
(0.01 < z < 0.07) early-type galaxies (ETGs) with prominent dust lanes, drawn from 
the Sloan Digital Sky Survey (SDSS). Two-thfrds of these 'dusty ETGs' (D-ETGs) 
are morphologically disturbed, which suggests a merger origin, making these galaxies 
ideal test beds for studying the merger process at low redshift. The D-ETGs preferen- 
tially reside in lower density environments, compared to a control sample drawn from 
the general ETG population. Around 80 per cent of D-ETGs inhabit the field (com- 
pared to 60 per cent of the control ETGs) and less than 2 per cent inhabit clusters 
(compared to 10 per cent of the control ETGs). Compared to their control-sample 
counterparts, D-ETGs exhibit bluer ultraviolet-optical colours (indicating enhanced 
levels of star formation) and an active galactic nucleus fraction that is more than an 
order of magnitude greater (indicating a strikingly higher incidence of nuclear activ- 
ity) . The mass of clumpy dust residing in large-scale dust features is estimated, using 
the SDSS r-band images, to be in the range lO^'^-lO^'^M©. A comparison to the total 
(clumpy -I- diffuse) dust masses - calculated using the far-infrared fluxes of 15 per 
cent of the D-ETGs that are detected by the Infrared Astronomical Satellite (IRAS) 
- indicates that only 20 per cent of the dust is typically contained in these large-scale 
dust features. The dust masses are several times larger than the maximum value ex- 
pected from stellar mass loss, ruling out an internal origin. The dust content shows no 
correlation with the blue luminosity, indicating that it is not related to a galactic scale 
cooling flow. Furthermore, no correlation is found with the age of the recent starburst, 
suggesting that the dust is accreted directly in the merger rather than being produced 
in situ by the triggered star formation. Using molecular gas-to-dust ratios of ETGs 
in the literature, we estimate that the median current molecular gas fraction in the 
IR^S'-detected ETGs is ~1.3 per cent. Adopting reasonable values for gas depletion 
time-scales and starburst ages, the median initial gas fraction in these D-ETGs is ^4 
per cent. Recent work has suggested that the merger activity in nearby ETGs largely 
involves minor mergers (dry ETG -I- gas-rich dwarf), with mass ratios between 1:10 
and 1:4. If the ZR^S'-detected D-ETGs have formed via this channel, then the original 
gas fractions of the accreted satellites are between 20 and 44 per cent. 

Key v^rords: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: interactions - galaxies: ISM - galaxies: peculiar. 
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1 INTRODUCTION 

As endpoints of the hierarchical build-up of mass, early- 
type galaxies (ETGs) possess a detailed fossil record 
of their progenitors. The stellar populations of ETGs 
encode the mass assembly history of galaxies over the 
lifetime of the Universe, and it is important that we 
gain a comprehensive understanding of their forma- 
tion and evolution. Whi le the ir red optical c olours 
' e.g. iBower. Lucev fc Eii3 Il992l: iBernardi e\~^ l2003l : 



iBell et al.l |2004| : | Faber et al 
enhancements (e.g 



2007 ) , high alpha- el ement 

. iThomas. Greggio fc Bendeil 1 19991 : 

iTrager et al.l l2000al lb[) an d obedience of a tight Fun- 



dame nt al Plane (e.g. Jorgensen. Franx fc KiaergaardI 



19961; Ivan Dokkum fc Franxl Il996l : ISaglia et al.l Il997l : 



Forbes. Ponman fc Brownl 1 19981 ) indicate that the bulk of 
their stellar populations (> 80 per cent) are old, the large 
scatter in the star-formation-sensitive ultraviolet (UV) 
colours of ETGs since z ~ 1 is interpreted as evidence 
for continuous low- l evel s t ar formation at late epochs 
llKavirai et all |2007|. l2008l: iMartin. O'ConneU fc HibbardI 



20091. 



also iFukugita et all 120041: lYi et al.l l2005l: 



20T1). 



Salim fc Rich! 
A coincidence between blue 



Jeong et al. 200 



20101 : ICrockett et"al] 
UV colours and mor- 



phological distu rbances indicates that t his star formation is 
merger driven fech weizer et al. 1990l: ISchweizer fc Seitzen 



ll992l : lKavira! 2010: Kaviraj et al.ll201ll 'l. Further more, the 
pauc i ty of major mergers at late epochs (see e.g . ISteward 
l2008l : |Jogee et al.ll2009l : iLopez-Saniuan et al.ll201lh strongly 
suggests that the sta r form ation is driven by minor mergers 
l|Kavirai et al.l |2009| . 120111 ) . Our current understanding of 
ETG evolution therefore indicates that their underlying 
stellar populations are old, having been rapidly built at 
high redshift (z > 1), while the evolution at late epochs 
{z < 1) is dominated by repeated minor-merger events that 
contribute < 20 per cent of their stellar mass at the present 
day. 

Notwithstanding the classical notion of ETGs be- 
ing dry, passively evolving systems, an extensive lit- 
erature has developed over the past few decades on 
the interstellar medium (ISM) in these systems. The 
dust and gas contents of very nearby ETGs have 

[Tubb3 1 19801: 



been studied 



.% 



several 
r 



authors (e.g . 



1985": 'Bcrtola! 



Hawarden et al. 1981; Sadler fc Gerhard 

1987; Ebnctcr, Davis & Djorgovski 1988; Knapp ct al. 1989; 



Goudfrooii fc de Jond 1 19951: Ivan Dokkum fc Franxl 1 19951: 



Knapp fc RupenI Il996l; iFaber et all Il997l: iTomita et all 
20od; iTran et al.1 I2OO1I: ICombes Young fc Bureaul l2007l : 



Calura. Pipino fc Matteuccill2008l : ivong et al.ll201ll ). While 
the galaxy samples and methodologies in these studies are 
varied, it is clear that the majority of ETG s in the very 
nearby Universe show evidence o f dust (e.g. iKnapp et al.l 
1 19891 : Ivan Dokkum fc Franxl Il995l ) , the dust masses being 
generally inconsistent with scenarios in which the dust is 
supplied purely by stellar mass loss (e.g. iMerluzzi 1998h . 
Coupled with the frequently observed morphological distur- 
bances and kinematical misalignments between gas and stars 
in these system s (e.g. Zcilingcr, Bcrtola & Gallctta 199S; 
Sage fc Galiettal 1 19931 : lAnnibah et al.1 120101 : iDavis et all 



201l|), it seems likely that a significant fraction of the in- 



terstellar matter has an external origin. 

The accumulating evidence for widespread minor- 



merger-driven star formation (e.g. lKavirai et al.ll2009l . 120111 ) 
makes the study of interstellar matter in ETGs all the more 
compelling. While the properties of the star formation have 
been quantified with a reasonable degree of precision, less 
is known about the fuel that drives this star formation. A 
study of the ISM of ETGs in modern observational surveys is 
therefore very desirable. In this paper we present a study of 
dusty ETGs (D-ETGs), drawn from Data Release 7 (DR7) o f 
the Sloan Digital Sky Survey (SDSS: iAbazaiian et al.ll2009l ). 
that exhibit prominent, extended dust features in their op- 
tical images. In addition to exploring the properties of the 
dust, the homogeneous nature of the spectrophotometric 
data from the SDSS allows us to systematically study the 
star formation, active galactic nucleus (AGN) activity and 
local environment of the D-ETGs, compare these systems to 
a statistically meaningful 'control' sample drawn from the 
general ETG population and put some of the results in the 
literature on a firmer statistical footing. We are also able to 
explore the properties of the gas in the minor mergers that 
are responsible for star formation in nearby ETGs, and pro- 
vide a more complete picture of the star formation activity 
in massive ETGs at late epochs. 

The plan for this paper is as follows. Section [2]describes 
the selection of D-ETGs from the SDSS and the construction 
of a control sample to which their properties are compared. 
In Section[3]we discuss the high incidence of disturbed mor- 
phologies observed in the D-ETG sample. In Section!?] we 
explore the local environment of these galaxies, and in Sec- 
tion [5] we employ UV-optical spectrophotometry to study 
their star formation and AGN activity. We explore the prop- 
erties and origin of the dust in Section [S] and conclude by 
summarising our findings in Section [7| 



2 DATA 

2.1 Sample selection and basic properties 

Galaxy Zoo (GZ: iLintott et aLllJOOJ . I2OIII ) is a citizen sci- 
ence project that has used more than 250 000 members of the 
general public to morphologically classify ~900 000 galaxies 
in the SDSS spectroscopic galaxy sample. The first incar- 
nation of this project [Galaxy Zoo 1 (GZl)] classified ob- 
jects in the SDSS galaxy population into broad morpholog- 
ical classes: early types, spirals and mergers. Galaxy Zoo 2 
(GZ2), an extension to GZl on which this study is based, 
has explored finer details of galaxy morphologies, such as 
morphological peculiarities (dust lanes, tidal disturbances, 
etc.). Note that the GZ objects are taken from the SDSS 
galaxy catalogue, which does not contain any quasars or 
type I AGN. Such broad-line objects are identified by a 
pipeline algorithm which compares the observed SDSS spec- 
trum of individual galaxies to the composite spectrum of 
IVanden Berk et al.l (|200ll l^l and do not appear in the galaxy 
catalogue. 

Around 19 000 GZ2 galaxies in the redshift range 0.01 < 
z < 0.1 are flagged as containing a dust feature by at least 
one GZ2 user. Each galaxy in this sample was visually re- 
inspected by S. Kaviraj and Y. S. Ting to determine (1) 



^ http://www.sdss.org/dr7/algorithms/redshift_type.html 
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Figure 1. Examples of ETGs with prominent dust features from the sample studied in this paper. Widespread morphological 
disturbances (shells and tidal tails) are visible in the galaxies. A colour version of this figure is available in the online version of the 
journal. Note that the faint morphological disturbances may be more readily visible on screen than in printed images. 



whether the galaxy really has a dust feature, (2) the mor- 
phology of the parent galaxy (ETG or non-ETG) and (3) 
whether the galaxy is morphologically disturbed. After this 
second layer of visual inspection we decided to restrict this 
study to the redshift range z < 0.07, because the unambigu- 
ous detection of a dust feature appears to become difficult 
beyond this redshift. Note that the identification of D-ETGs 
is not expected to be dependent on galaxy luminosity. As 
long as the dust is backlit by the galaxy, it is readily visi- 
ble. This process yields a final sample of reliably classified 
352 D-ETGs, which are the subject of this paper. They rep- 
resent 4 per cent of the ETG population (restricted to the 
same redshift and magnitude ranges as the D-ETG sample). 
Fig. [l] shows several typical examples of our D-ETGs. 
Fig- shows the basic properties of the D-ETG sample: the 
redshift, absolute r-band magnitude and stellar mass distri- 
bu tions. Stell a r mas ses are calculated using the calibrations 
of iBell et al.l l|2003l ) . All magnitudes shown in this paper 
have been if-correctcd, using the latest versio n of the pub- 
lic KC ORRECT code (,Blanton ct al., ,2003. : .Blanton fc Roweia 
l2007l ). 



2.2 Control sample 

An important objective of this paper is a systematic com- 
parison between the properties of D-ETGs and those of a 
control sample of galaxies from the general ETG population. 
We create a control sample of ~4000 ETGs, using galaxies 



from GZ2 which have an early-type 'vote fraction' greater 
than 0.7, i.e. more than 70 per cent of the users that classi- 
fied each of these objects fiagged it as an ETG. The control 
sample is selected to be 'dustless', i.e. no GZ user fiagged 
any of these galaxies as having a dust lane. The control 
sample is then visually inspected to remove late-type con- 
taminants (the contamination rate is very low, around 3 per 
cent) and flag control ETGs that have morphological distur- 
bances. The control sample is further constructed so that it 
has the same redshift and magnitude distributions as the 
D-ETG population shown in Fig. [21 



3 DISTURBED MORPHOLOGIES: A MERGER 
ORIGIN FOR DUSTY ETGS 

The visual inspection described above indicates that ~65 
per cent of D-ETGs show clear morphological disturbances, 
at the depth of the standard SDSS images. In comparison, 
the control sample (described in Section 12. 2p shows a dis- 
turbed fraction of only around 6 per cent. The fraction of 
disturbed objects in the D-ETG sample is therefore strik- 
ingly high, indicating that D-ETGs are the products of re- 
cent interaction events and, in agreement with the literature 
(e.g. Goudfrooij ct al. 1994), that the formation of the dust 
is likely to be the result of these interactions. Indeed the 
D-ETGs in this sample are likely to represent objects like 
Centaurus A in the making and are thus ideal test beds for 
studying the merger process in the low-redshift Universe. 
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Figure 3. Cumulative distribution of host DM halo masses for 
D-ETGs and the ETG control sample. Compared to the con- 
trol sample, D-ETGs preferentially inhabit lower density envi- 
ronments. Almost 80 per cent of the D-ETGs reside in the field 
(compared to 60 per cent of the control ETG population) and less 
than 2 per cent of the D-ETGs reside in clusters (compared to 10 
per cent of their control-sample counterparts). A colour version 
of this figure is available in the online version of the journal. 



isolated galaxies. This catalogue provides estimates of the 
masses of the host dark matter (DM) haloes of individual 
SDSS galaxies, which are related to the traditional classifica- 
tions of environment ('field', 'group' and 'cluster'). Cluster- 
sized haloes typically have masses greater than lO^^M©, 
while group-sized haloes have masses between 10^^ and 
10^'' Mq. Smaller DM haloes constitute what is commonly 
termed the field. 

In Fig. |3]we present the cumulative distribution of host 
DM halo mass for both the D-ETGs and the ETG con- 
trol sample. We find that, compared to the control ETG 
population, D-ETGs preferentially inhabit lower density en- 
vironments. Almost 80 per cent of D-ETGs reside in the 
field (compared to 60 per cent of the control ETGs) and less 
than 2 per cent reside in clusters (compared to 10 per cent of 
their control-sample counterparts). This is consistent with 
the idea that the D-ETGs are merger remnants, since the 
high rela tive velocities in cluster env ironments make mergers 
unlikely l|Binnev fc T^emainelll987^ . 



Figure 2. Distributions of redshift (top), r-band absolute mag- 
nitude (middle) and stellar mass (bottom) of the D-ETG sample 
in this study. As described in Section l2.2l the control samples to 
which these galaxies are compared are matched to the shape of 
the distributions presented here. 



4 LOCAL ENVIRONMENT 

We begin by probing the local environment of our D- 
ETGs, by cross-matching this sam ple with the SDS S en- 
vironment catalogue constru c ted b y lYang et al.l (|2007l ) and 
lYang. Mo. fc van den BoschI ()2008l ). who use a halo-based 
group finder to separate the SDSS into over 300 000 struc- 
tures with a broad dynamic range, from rich clusters to 



5 STAR FORMATION AND AGN ACTIVITY 

Fig. |4]shows the [g — r) and (NUV— r) colour distributions 
of the D-ETGs compared to those of their control-sample 
counterparts. The optical g and r magnitudes are taken 
from the SDSS _York_et aL (;20.00,). The NUV magnitudes 
are from the Galaxy Evoluti on Explorer (GALEX) space 
telescope jMartin et al.ll2005h. The GALE X NUV filter is 
centred at ~2300 A JMorrissev et al■ll2007^ . We use the rec- 
ommended total magnitudes in the AB systems for both 
the SDSS and GALE X photometry (|Stoughton et al.ll2002l : 
iMorrissev et al.ll2007l '). 

We find that the D-ETGs are redder than their control 
counterparts in the optical {g — r) colour, which is consis- 
tent with the visibly high levels of dust present in these 
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Figure 4. Top: optical {g — r) colour distributions of the D-ETGs 
and the control ETGs. Unsurprisingly, the D-ETGs are found to 
be redder than their control counterparts, due to the presence of 
visibly large amounts of dust. Bottom: (NUV— r) colour distri- 
butions of the D-ETGs and the control ETGs. The D-ETGs are 
bluer in the UV-optical colour, presumably due to recent star for- 
mation (see text for details). Median values are shown using the 
dashed vertical lines. A colour version of this figure is available 
in the online version of the journal. 



systems. Nevertheless, the D-ETG population is bluer than 
the control sample in the (NUV— r) colour, with a median 
offset of ~0.5 mag. Given that the UV wavelengths are sev- 
eral times more sensitive to young stars than the optical 
wavelengths (see e.g. fig. 1 in lKavirai et al.ll2009l ). the bluer 
UV-optical colours are likely to be driven by recent or ongo- 
ing star formation in the D-ETG population. It is important 
to note that the large-scale structure of the dust in the D- 
ETGs is clumpy and does not entirely obscure the galaxy 
core (see Fig. [ij, where the young stars (w hich dominate 
the UV flux) are typically concentrated (e.g. Ferreras et al.l 



I2OO9I : iJeong et al.ll2009l : iPeirani et ai]|201l]| ; ISuh et al.ll2o7of )^ 
We proceed by probing the dominant gas ioniza- 
tion mechanism in the D-ETG sample. The presence 
of dust implies the availability of gas which, in turn, 
makes the presence of star formation and nuclear activ- 
ity in these galaxies likely. We use the standard optical 
emission- line ratio analysis of iKewlev et al.l (|2006l . see also 



Table 1 . Comparison of star formation and AGN activity, derived 
using optical emission-line ratio diagnostics, in the D-ETGs and 
their control-sample counterparts. 



D-ETGs 
(per cent) 



Control ETGs 

(per cent) 



/(D-ETG)//(Ctrl) 



Star forming 


7.2 


2.1 


Composite 


22.8 


2.8 


Seyfert 


19.9 


1.7 


LINER 


37.5 


10.4 


Quiescent 


12.8 


83.0 



3.4 
8.1 
11.7 
3.6 
0.2 



Baldwin. Phillips fc TerlevichI Il98ll : IVeilleux fc Osterbrockl 
19871 : iKauffmann et al.l |2003| ) which exploits the [Nil] /Ha 
and [Oin]/H/3 ratios to probe the presence of type II AGN 
in our D-ETG sample. The optical emission-line fluxes are 
calculated from the SDSS spectru m of each object using 
the publicly available code GANDALF (|Sarzi et al.l|2006l ). Note 
that the SDSS spectra are measured in 3-arcsec fibres (^ 4 
kpc at z = 0.035, the median redshift of our D-ETG sam- 
ple) centred on the galaxy. Galaxies in which all four emis- 
sion lines are detected with a signal-to-noise ratio greater 
than 3 are classified as either 'star forming', 'composite' (i.e. 
hosting both star formation and AGN activity), 'Seyfert' 
or 'LINER'. Galaxies without a detection in all four lines 
are classified as 'quiescent'. The fraction of galaxies in each 
emission-line class is listed in Table [iJ for both the D-ETG 
and control samples. 

We find that the fractions of D-ETGs classified as star 
forming and composite are, respectively, factors of 3 and 8 
larger than that for the control ETGs. This suggests en- 
hanced levels of star formation in the D-ETG population 
compared to its control counterpart, consistent with the 
UV-optical colour analysis above. The fraction of LINERs 
in the D-ETGs is also greater by a factor of 3, plausibly 
due to LINER-like emission from merger-driven shocks or 
evolved stellar populations (such as post-asymptotic giant 
bran ch stars) that ar e associated with the star formation 
(e.g. lSarzi et al.ll2010l ). Table [l] indicates that, in addition to 
enhanced star formation and LINER emission, the D-ETGs 
exhibit a strikingly higher Seyfert fraction, more than an or- 
der of magnitude larger than in their control-sample counter- 
parts. The fraction of quiescent galaxies is correspondingly 
an order of magnitude smaller. Therefore, the presence of 
prominent dust features in an ETG significantly increases 
its chances of hosting both more star formation and more 
nuclear activity, with nuclear activity being the dominant 
gas ionization mechanism in these systems. 



6 PROPERTIES OF THE DUST AND 
ASSOCIATED GAS 

6.1 Clumpy dust 

The dust in D-ETGs is comprised of both clumpy structures, 
which are easily visible in the optical images, and diffuse 
dust, which is intermixed with the stars and is essentially 
invisible in the optical imaging. We begin by calculating the 
dust contained within the clumpy features (M^^), follow- 
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ing a metho d similar to that employed by several previous 
studies (e.g. Sadler fc Gerhardlll985l : Ivan Dokkum fc Franxl 
1 19951 : iTran et al.ll200ll ). Thedumpv dust mass is given by 



Mf^ 



{Av)T.Tl 



(1) 



where {Ay) is the mean V-band extinction in magnitudes, 
E is the total area of the dust features and Ty is the mass 
absorption coefficient [j 

Following these previous studies, Av is measured as fol- 
lows. Ellipses are first fitted to the isophotes of the SDSS 
r-band images of each galaxy, using the IRAF task ELLIPSE, 
with the dust features masked out. An absorption map is 
then created, using the ratio of the flux in the observed im- 
age (-Fobs) to the flux in the model fit (-Fmod). The optical 
depth is r = — ln(Fobs/Fniod) and Ar = 1.09r l|Tran et al.l 
I2OOII ). We con vert Ar to Ay us ing the Galactic dust law: 
Ay = 1.33^r (|Tran et al.ll200ll '). Note that the quality of 
the ELLIPSE fits are unaffected by the faint tidal features 
that are found on the outskirts of the D-ETGs. 

Several studies have shown that the extinction curves 
in ETGs with du st lanes are similar to their Galactic coun- 
terpa r t (see e.g. Sadler fc Gerhard! 19851 : iGoudfrooii et al 
19941: ISahu. Pandev fc Kembhavil Il998l : iFinkelman et al 



20081 ). which indicates that the properties of the dust in 



these systems are similar to that in the Galaxy. Therefore, 
it appears reasonable to employ the Galactic value for the 
mass absorption coefficient (i.e. Fv = 6 x 10~^ niag kpc^ 



m; 



although we note that a more accurate calculation 



requires a measurement of Ty in our individual D-ETGs. 
The formal error on M^^ (driven, in this case, only by the 
error in {Ay)) is calculated using the standard error propa- 
gation formula. We note that the errors in other quantities 
that appear later in this study (e.g. M^^ and gi) are also 
calculated using the standard error formula. 

While the formal error on M^^ is only ~0.055 dex, this 
derivation assumes that the dust is placed in a simple fore- 
ground screen in front of the light source; however, in actual 
fact the dust is embedded in the galaxy. The use of this sim- 
plifying assumption induces an additional uncertainty o f up 
to ~0.3 dex in the dust mass calculation (e.g. Tomita et al.l 
I2OOOI : iTranet al.l 120011 ). The values derived here (see the 
top panel of Fig. [SJ are in the range IO^'^-IO'^'^Mq, in 
good agreement with the range of clumpy dust masses 
derived using similar techn i ques i n the literature (10^- 
IO'^M m : [Sadler fc GerhardI Il985l: IGoudfrooii fc de Jond 
I995I: Iraji Dokkum fc Franxl Il995l: iToniita et al.l I2OO0I : 



Tran et al.ll200ll : IWhitaker fc van Dokkumll20oi r 

In the middle panel of Fig. [S] we plot the clumpy 
dust mass of the D-ETGs against their total blue lumi- 
nosity [Lb)- Lb is calculated from the absolute (total) 
-B-band magnitude, after transforming from the absolute 
SDSS (to tal) fl-band magnitu de using the relation nis = 
mg -\- 0.3 JRodgers et al.ll2006l ). We use Lb to facilitate d i- 
rect comparison with previous studies (e.g. lMerluzzilll998h . 
The dashed line indicates the maximum expected contribu- 
tion to the dust content from stellar mass loss alone (e.g. 



■^ The right-hand side of equation Q is derived as follows. If 5 is 
the area of a single pixel and Ai is the extinction in pixel i, then: 



[Ml 



E = nS, M^L 



- SA2 + M3 

= (E/n)[Ai- 



- • • ■ -f <5A„]r- 

A2 + As + --- 



Since the total area 





10.0 10.5 

Log L, [Lg] 



^ 6 




-2 -1 

Log age of recent storburst [Gyrs] 



Figure 5. Top: distribution of clumpy dust mass (M^^) in the 
D-ETGs, derived from their optical r-band images (see text for 
details). Middle: clumpy dust mass versus blue luminosity of the 
D-ETGs. Bottom: clumpy dust mass versus age of the recent star- 
burst, estimated from UV-optical colours. Note that, while the 
formal error on M^^ is ~ 0.055 dex (driven by the uncertain- 
ties in Ay), the uncertainties in the assumed dust distribution 
induces an additional error of up to 0.3 dex (see text for details). 
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Faber fc Gallagheij 1 19761 : iKnapp. Gunn fc Wvnn-Williamsl 
1992f). assuming that the dust is depleted b y sputtering 
'e.K- Spitzer 1978; Drainc " " ' ' -«-« . -.-. 



lust IS depleted b y sputtering 
Sal peteilll979 a"b': Dwck'l98d; 



iDwek. Rephael i & Mather 1990,) over a m aximal destruc- 
tion time-scale of ~ 10^'^ yr (IFaber fc Ga llagher 1976). In 
agreement with previous studies fe.g. lGoudfrooii et al.iil994l : 
iMerluzzillTggi ). we find that (1) D-ETGs typically have 
clumpy dust masses that are well in excess of the maximum 
equilibrium mass expected from stellar mass loss alone and 
(2) there is no correlation between the optical dust masses 
and Lb ■ It is worth noting that Lb shows a correlation with 
the total X-ray lumi nosity (Lx) of ETGs, which traces their 
hot gas content (e. g. Caniza res. Fabbiano fc Trinchierl 19871 : 
iMathew fc BrighentL2003. : .Mulchaev fc Jeltemall2010l ). The 



lack of a correlation between the clumpy dust mass and Lb 
suggests that the clumpy dust is not produced through cool- 
ing from the hot gas rese r voirs t hat may be present in these 
ETG s (see also iForbed Il99ll : iBregman. Hogg fc Robert^ 
1 19921 ). Recall also that the D-ETGs studied in this paper 
do not typically inhabit clusters, so that cluster-scale cool- 
ing flows are not relevant to these systems. 

The dust properties derived above indicate that the 
clumpy dust content of the D-ETGs is not created by in- 
ternal processes such as stellar mass loss or gas cooling. 
Together with the high incidence of morphological distur- 
bances observed in the D-ETG sample (see Section [S]), this 
indicates that the clumpy dust observed in these galaxies is 
likely to be associated with the recent merger events that 
created these objects. 

Since star formation is clearly enhanced in D-ETGs, it 
is worth exploring whether the dust is produced in situ by 
this star formation or whether it is directly accreted in the 
mergers. In the bottom panel of Fig. [S] we plot the clumpy 
dust mass of D-ETGs against the age of the recent starburst 

ST), which has been estimated using the UV-optical colour 
Kaviraj|2010|): 



logT [Gyrs] 



0.6"-''^ X [(NUV - 



(9 



1.73=" 



Briefly, (see sections 4.2.1 and 4.3 in lKaviraill20ld . for de- 
tails), T is one of several star formation history (SFH) 
parameters derived for individual ETGs by fitting their 
GALEX /ST>SS UV-optical photometry to a library of syn- 
thetic photometry, generated using a large collection of 
model SFHs. The model SFHs assume two instantaneous 
bursts, the first fixed at z — 3 (since the bulk of the stars 
in ETGs are old), while the second (which represents the 
recent star formation episode) is allowed to vary in age 
(T) and mass fraction (realistic ranges of dust a nd m etal- 
licity are applied to the model SFHs). iKavirail (|2010|) ex- 
tract marginalized values for the free parameters including 
T and find a strong correlation between the marginalized 
values of T and the relatively dust-insensitive double colour 
[(NUV— m) — {g— z)]. This relation (equation [21) allows us to 
estimate the age of the recent starburst in ETGs from the 
observed value of [(NUV— ii) — (g — z)]. 

Fig. [S] shows that the clumpy dust mass does not cor- 
relate with T. The apparent decoupling between the age of 
the recent (merger-induced) star formation episode and the 
clumpy dust content in D-ETGs favours a scenario in which 
the dust is directly accreted during the merger and not cre- 
ated in situ by the triggered star formation. 




5.0 5.5 6.0 

Log K' [Me] 

Figure 6. Comparison of the clumpy dust masses {M^^), derived 
from the optical r-band images, with the dust masses derived from 
IRAS fluxes M^^^S T^e JRAS'-derived dust masses are typically 
a factor of 4.5 larger than their clumpy counterparts. Note that 
the error bar in M^^^ is the error from equation JS}- 



6.2 Diffuse dust and the total dust content 

Since the diffuse dust mass in the ISM cannot be estimated 
from the optical images, the clumpy dust masses calculated 
above are likely to be lower limits to the total dust masses. 
It is thus desirable to quantify the total (diffuse + clumpy) 
dust content of these systems, at least to gain an under- 
standing of how much dust is not contained in the clumpy 
features visible in the optical images. A useful route to es- 
timating the total dust content is to use the far-infrared 
spectrum. Around 15 per cent of our D-ETGs have infrared 
data from the Infrared Astronomical Satellite (IRAS), which 
we use to estimate the 'total' dust mass in these galaxies. 
It should be noted, however, that IRAS is insensitive to 
(2) very cold dust, with a temperature less than ~20 K (e.g. 
iHenkel fc Wikhnd|[l997l ) . Thus, while the IRAS fluxes are a 
much better estimator of the total dust mass than the opti- 
cal images, more accurate estimates require imaging by new 
instruments such as Herschel, that probe the far-infrared 

spectrum beyond 100 /^m. 

Following pa st studies fe.g. lGoudfrooii fc de J ong 199J; 
iTranet al.ll200ll ). the IRAS-derweA dust i nasses (Mj^^^) 

' tj Il982l : 



are calculated using the formula (see e.g. Schwart 



iHildebrandl 19831 : lRowan-Robinsonlll986l : iKwan fc Xij|l992l ) 



M. 



IRAS 



= 5.1 X lQ-''S^D^Xt[e 



2\4r (1.44xlO''/A„Ta) 



1] Mq, (3) 



where the dust temperature (Td) is estimated using the 
IRAS flux ratio Seo/Swo and an emissivity law that varies 
as A"'^, 5*1, is the IRAS flux density in mjy at wavelength 
i^, D is distance to the galaxy in Mpc and \v is the wave- 
length in microns. We use v = 100 /xm for our calculations. 
Fig. m shows the IRAS-AexiveA dust masses plotted against 
the clumpy dust masses, for the subset of D-ETGs that 
have IRAS detections at 60- and 100-^m. The IRAS fluxes 
are taken from the Imperial IRA S-FSC Redshift Catalogue 
ijWang fc Rowan- Robinson! 120091 ). The median discrepancy 
between the clumpy dust mass and the IRAS dust mass is a 
factor of 4.5. In other words, only ~20 per cent of the dust 
is typically contained in clumpy features. Nevertheless, the 
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Figure 7. Derived current (top) and initial (bottom) molecular 
gas fractions (i.e. molecular gas mass divided by the stellar mass) 
in the subset of D-ETGs that are detected by IRAS (see text for 
details). Median values of the distributions are indicated using 
solid (red) lines. 



IRAS dust masses show the same lack of correlation with 
blue luminosity and starburst age as the clumpy dust masses 
do, leaving our conclusions above unchanged. 



6.3 Cold molecular gas 

Give n the close co rrespondence between dust and cold gas 
(e.g. lKnapdll999l ). it is instructive to explore the proper- 
ties of the cold gas reservoirs of our D-ETGs. We choose 
to focus on molecular gas (i.e. H2, derived from CO ob- 
servations) because of its relevance to star formation (e.g . 



Wong fc B litz"2002l: ISchavell2004l : iKrumholz fc McKed[200g : 
Lerov et a l. 2008) and the better availability of molecular 
gas data on ETGs in the literature compared to its atomic 
counterpart (Hi). We note first that previous studies have 
typically assumed that molecular gas-to-dust ratios (G/D 
values) in ETG s arc similar to that i n the Galaxy (~150 ; 
see e.g. lSpitzeij|T97 8: Hildcbrand 198J: lDraine fc LedflQSi ). 
However, if the dust is directly accreted during recent merg- 
ers (c.f. Section fG.ip . then the Galactic G/D is inappropriate 
for all ETGs, unless the accreted companions are all Milky 
Way like. In fact, recent work (e.g. lKavirai et al]|2009l . l201ll : 



ICrockett et al.ll2011h has shown that a large proportion of 
the merger activity experienced by nearby ETGs involves 
interactions with small gas-rich companions (rather than 
Milky Way like systems), implying that the Galactic G/D 
value may not be appropriate for our D-ETG sample. 

To obtain a more realistic conversion between dust 
and molecular gas for our D-ETGs, we appeal to mea- 
surements of molecular G/D values in ETGs that are 
available in the literat u re. For 76 ETGs st u died by 
IWiklind fc Henke| (Il989l ): iThronson et all (Il989l): iGordonI 
('1991'): 'Logs ot al.^ ('l99l'):'Wang. Kennev fc Ishizukil (|l992f ) 
and Wiklind fc Hcnkcl ( 1995), the average molecular G/D is 
1190 ±410. Thus, the molecular G/D values appropriate for 
our D-ETGs are likely to be almost an order of magnitude 
higher than the value in the Galax;y. We note here that the 
measured G/D values in ETGs are indeed around t he values 
typical of late-type dwarfs (e.g. lYoung et al.l ll98!l). further 
supporting the view that the interactions that produce the 
D-ETGs probably involve the accretion of gas-rich satellites. 
In the top panel of Fig. [T] we use the IRAS dust masses 
to estimate the current molecular gas contents, assuming a 
G/D of 1190 ± 410. The median molecular gas fraction (i.e. 
molecular gas mass divided by stellar mass) is ~1.3 per cent. 
The error bars incorporate both the uncertainty on the G/D 
value and the errors in the IRAS-deiived dust masses. It is 
worth noting that this value is well within the distribution 
of molecular gas fractions found bv I Yong et al.l (2011 ) in a 
volume-limited sample of normal ETGs in the ATLAS'^'^ 
survey of very nearby ETGs. 

Since star formation is ongoing in these galaxies, we ex- 
pect the current molecular gas masses calculated above to 
be lower than the initial molecular gas contents. Given a 
characteristic exponential time-scale (r^) for the gas deple- 
tion (which is also the star formation time-scale) , the current 
molecular gas fraction (go) is related to the initial value (gt) 
by 



50 = pie 



(-T/t^) 



(4) 



where T is the look-back time to the onset of star forma- 
tion (i.e. the starburst age). Using equation ((2| to estimate 
r, and adopting a value of 0.1 Gyr for the characteristic 
time -scale of the recent st ar formation episodes in ETGs 
(see ISchawinski et al.l 120071 ) , we estimate gt in each of our 
D-ETGs. These initial molecular gas fractions are presented 
in the bottom panel of Fig. [T] The estimated initial gas frac- 
tions are small, wit h a median value of around 4 per cent. Re- 
cent work (see e.g. iKavirai et al.ll2009l . I2OIII ) has suggested 
that the merger activity (and associated star formation) in 
nearby ETGs largely involves minor mergers between dry 
ETGs and gas-rich satellites with mass ratios between 1:4 
and 1:10. If the IRAS-detected D-ETGs indeed formed via 
this channel, then a typical initial gas fraction of 4 per cent 
implies that the gas fractions of the accreted satellites are 
likely to be between ^^^20 and ~44 per cent (assuming that 
the ETGs themselves are gas free). It is worth noting that 
these empirical results are consistent with recent numerical 
work which suggests that the satellites fuelling the star for- 
mation in nearby ETGs should have molecul ar gas fractions 
greater than 20 per cent (|Kavirai et al.ll2009l ). It should also 
be noted that the molecular gas analysis is only represen- 
tative of the 15 per cent of D-ETGs that are detected by 
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IRAS. A more robust statistical analysis requires infrared 
data on the entirety of the D-ETG sample studied here. 



7 SUMMARY 

We have explored the properties of dust and associated 
molecular gas in ETGs, by studying 352 nearby (0.01 < z < 
0.07) ETGs with prominent dust features in the SDSS DR7. 
The galaxies were identified through direct visual inspec- 
tion of SDSS images using the GZ2 project. Our analysis 
has focused on comparing the local environment, star for- 
mation and nuclear activity of these D-ETGs to those of a 
control sample drawn from the general ETG population. We 
have then explored the properties and origin of the dust and 
associated molecular gas in these systems, and studied the 
processes that lead to their formation. 

Two-thirds of the D-ETGs show morphological distur- 
bances at the depth of the standard SDSS images, suggesting 
a merger origin for these galaxies. D-ETGs reside preferen- 
tially in low-density environments, with 80 per cent inhabit- 
ing the field (compared to 60 per cent of their control-sample 
counterparts) and less than 2 per cent inhabiting clusters 
(compared to 10 per cent of their control-sample counter- 
parts). D-ETGs show bluer UV-optical colours compared 
to the control ETG population, indicating the presence of 
enhanced levels of star formation activity. An analysis us- 
ing optical emission-line ratios indicates that the fraction of 
objects classified as star forming and composite are, respec- 
tively, factors of 3 and 8 higher in the D-ETGs compared to 
their control-sample counterparts. The fraction of Seyferts is 
more than an order of magnitude higher, while the fraction 
of galaxies classified as quiescent is correspondingly an order 
of magnitude smaller in the D-ETG population compared 
to that in the control sample. The presence of prominent 
dust features in an ETG therefore significantly increases its 
chances of hosting both star formation and nuclear activ- 
ity, with nuclear activity being the dominant gas ionization 
mechanism. 

We have explored the dust content of our D-ETGs, 
both in terms of the clumpy dust mass contained in the 
large-scale dust features (visible in the optical SDSS im- 
ages) and in terms of the diffuse dust that may exist in 
the ISM. The clumpy dust masses span the range 10*'''- 
10^'^Mq and are typically several factors higher than the 
maximum dust contribution expected from stellar mass loss 
from the old, underlying stellar populations that dominate 
ETGs. The clumpy dust mass does not correlate with the 
blue luminosity of the galaxies, indicating that galactic-scale 
cooling flows are unlikely to be responsible for the dust con- 
tent. Similarly, no correlation is found between the clumpy 
dust content and the age of the recent starburst, suggesting 
that the dust is not produced in situ from the star forma- 
tion triggered by the merger event, but is rather accreted 
directly during the merger. 

The clumpy dust masses are strictly lower limits to the 
total dust masses because they do not include diffuse dust 
which may be intermixed with the stars and is essentially 
invisible in the optical images. In the 15 per cent of our 
D-ETGs that are detected by IRAS, we have estimated the 
'total' (clumpy + diffuse) dust mass using the IRAS 60- and 
100- /im fluxes. The average discrepancy between the clumpy 



and /7?j4S'-derived dust mass in these D-ETGs is a factor of 
4.5. In other words, only ~20 per cent of the dust mass in 
these objects is contained in the large-scale dust features 
that are observed in the optical images. Nevertheless, the 
IRAS dust masses show the same lack of correlation with 
galaxy properties and starburst ages as their clumpy coun- 
terparts. 

The correspondence between dust and cold gas allows 
us to explore the properties of the cold gas reservoirs in our 
D-ETGs. Restricting ourselves to the IRAS -detected galax- 
ies (in which we have a reasonable estimates of the total 
dust mass) , we employ the average molecular G/D of ETGs 
in the literature ('^ 1190) to estimate the current molecular 
gas fractions of these objects. The median value of the cur- 
rent gas fraction is ~ 1.3 per cent. Given the current molecu- 
lar gas fraction, an estimate of the age of the starburst and 
a characteristic gas depletion (star formation) time-scale, 
we have estimated the initial molecular gas fraction in each 
IRAS-detected D-ETG, with a median value of ^4 per cent. 
Given recent work that suggests that the merger activity in 
nearby ETGs largely involves minor mergers with mass ra- 
tios between 1:10 and 1:4 (e.g. iKayiraj et al. 2009), a typical 
initial molecular gas fraction of 4 per cent suggests that the 
gas-rich satellites accreted by the IRAS-detected ETGs are 
likely to have had molecular gas fractions between ~20 and 
44 per cent (assuming that the ETGs themselves are gas 
and dust free). 

While recent UV-optical studies have measured the con- 
tinuous (minor-merger-driven) star formation in ETGs over 
the last 8 billion years, less is known about the dust and 
associated gas that drives this star formation. This paper 
has offered insights into the ISM of nearby ETGs that are 
merger remnants and, in particular, the initial gas masses 
plausibly brought in by the satellites that drive the star for- 
mation. The galaxy sample studied here provides an ideal 
basis for future studies of D-ETGs using new instruments 
such as Herschel and the Atacama Large Millimeter Array 
(ALMA) to probe the properties of the ISM of nearby ETGs 
and achieve a fuller understanding of the gas that drives the 
star formation in these galaxies at late epochs. 
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